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We have investigated matter enhanced neutrino oscillations with a mantle-core-mantle 
step function and a realistic Earth matter density profile in both a two and a three 
neutrino scenario. We found that the realistic Earth matter density profile can be well 
approximated with the mantle-core-mantle step function and that there could be an 
influence on the oscillation channel — » Pt due to resonant enhancement of one of the 
mixing angles. 

1. Introduction 

Step function approximations o£lii£LJilarth matter density profile have been used by 
several authors, see e.g. Refs. tralj'Q, for the analysis of matter effects in neutrino 
oscillation experiments. In order to investigate the quality of this approximation, 
we have developed methods to numerically calculate the transition probability am- 
plitudes in a realistic matter density profile. Calculations were performed in a two 
and a three neutrino framework. The neutrino traveling path length L was divided 
up into several small intervals. In each of these intervals, in which we assumed 
constant matter density, we computed the effective mixing parameters. Combining 
the probability amplitudes of each interval leads to the total probability amplitude 
from which we obtained the transition probabilities. 

We have used three different matter density profiles, two model-like and one 
realistic: a constant profile just to show the resonance behavior, a step function 
describing the mantle-core-mantle structure of the Earth, and a realistic Earth 
density profile i. See Fig. |l| for the different density profiles. For the mantle-core- 
mantle density profile, we chose p = 5 g/cm.^ for the mantle, p — 12 g/cm.^ for the 
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core, and a core-width equal to half of the Earth's diameter. 

In our study we used the mass squared difference as suggested by atmospheric 
neutrino data from the Super-Kamiokande Collaboration obtained within a two 
neutrino analysis B 

Am^ ~ 3.2 • 10"^ eV^ 

For the constant density profiles in the case of two neutrino oscillations, we con- 
sidered the well-known Mikheyev-Smirnov-Wolfenstein (MSW) Q resonance condi- 
tion 

cos2^-— ^=0, A^A{p,E,) ^ ^^^""^ pE,, (1) 
Am^ rriN 

where 9 and Am^ are the usual mixing parameters for two neutrino oscillations, Gp 
is the Fermi weak coupling constant, Y (~ 1/2) is the average number of electrons 
per nucleon, m^r is the nucleon mass, p is the matter density, and is the neutrino 
energy. Using Eq. (|l|) when approaches zero, we obtained the resonant energies 
E,y ~ 3.6 GeV for the core and E^, ~ 8.5 GeV for the mantle. 
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Fig. 1. The matter density profiles. Realistic matter density profile (solid curve) and mantle-core- 
mantlc step function (dashed curve). 



2. Two Neutrino Calculation 

As mentioned above, we used Ato^ = 3.2 • 10^^ eV^ for the mass squared difference. 
Being interested in oscillation enhancement, we chose a small mixing angle 9 = 0.1. 
Using these mixing parameters, we calculated the transition probability P(i^e i^x), 
where x could be p, or r, for neutrinos going through the whole Earth (nadir angle 
zero). The results for the different matter density profiles discussed above are 
displayed in Fig. |^ as a function of the neutrino energy Ey. 

In the absence of matter (Fig. |2| (a)), there is a constant maximal probability 
of sin^ 26* without any enhancement. For constant matter densities (Fig. || (b)). 
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Fig. 2. The transition probability P(i/e — > u^:) as a function of the neutrino energy i?^. (a) zero 
matter density, (b) constant matter densities (mantle [solid curve] and core [dashed curve]), (c) 
mantle-core-mantle step function, and (d) realistic Earth matter density profile. 



mantle and core, we obtained oscillation enhancements at the expected neutrino 
energies of ~ 8.4 GeV and ~ 3.5 GeV, respectively. Figure ^ (c) shows the 
spectrum for the mantle-core-mantle model. The resonance spectrum is interpreted 
in the following way i'i'i: The main resonance is due to interference and is therefore 
not ]VISW-like. The shoulder to the left on top of the main peak is caused by the 
core and the peak to the right is caused by the mantle. Both these resonances 
are typically MSW-like. As seen from Figs. || (c)-(d), there is no large qualitative 
difference between the matter enhancements obtained using the step function and 
the realistic Earth density profiles. It turns out that the results for the step function 
profile strongly depend on the widths of the mantle and the core. 

The two neutrino model is not suitable for studying matter effects on the —>■ 
channel. To obtain more general results, we performed this investigation also in a 
three neutrino framework. 

3. Three Neutrino Calculation 

We used a hierarchal mass scheme Am^i <C Amgg ~ Ar7i§]^, where the large 
mass squared difference should describe atmospheric oscillations and the small mass 
squared difference should be responsible for solar oscillations. As in the two neutrino 
analysis, we again chose the large mass squared difference as Amgj = 3.2 • 10^"^ cV^. 
The small mass squared difference Am2i was set to zero under the assumption that 
these transitions play no role at length scales of the diameter of the Earth and 
smaller. 

We also used the standard parameterization for the neutrino mixing matrix with 
the mixing angles 6*12, ^13, and 623- To accommodate the atmospheric neutrino 
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result, we chose ^23 = 7r/4 giving the following vacuum transition probabilities 





I'e) 


= 1 sin2 2^?i3sin^^"J^^, 
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The parameters relevant for oscillations involving Am32 are ^13 and 023- The 
magnitude of 6*12 plays no crucial role even in presence of matter effects. For the 
numerical calculations we set ^12 — but also a maximal angle would not change 
our obtained results. 

Reactor experiments, the solar neutrino deficit and the atmospheric neutrino 
anomaly all indicate that the mixing angle ^13 should be small. The most stringent 
upper bound of sin^ 2^13 = 0.10 (valid for Am^ > 0.7 • 10""^ eV^) is coming from 
the CHOOZ experiment i. As in the two neutrino calculation we chose 6*13 = 0.1, 
which is obviously below the CHOOZ upper bound. 

The oscillation probability P{i'e ~> i^e) is given by an effective two neutrino 
formula with ^13 being the relevant mixing angle. Thus, the two neutrino model 
discussed above should give a good approximation for the oscillation enhancement 
obtained in this channel. The resonant energies in the three neutrino model are 
expected to be approximately the same as calculated in the two neutrino model. 
This can be seen in Fig. || (a), where the oscillation probability 1 — P{ve ^e) is 
shown as a function of the neutrino energy E^ for the realistic Earth density profile. 
Figure || (a) is very similar to the corresponding two neutrino case, see Fig. I (d). 

From Eq. (|^) it can be seen, that in the non-resonant region, electron neutrino 
disappearance is caused in equal amounts by the transitions Ve — t'^ and ^ Vt- 
This remains true also in the presence of matter as is shown in Fig. ^ (c). 

Due to the resonance enhancement of ^13, we expect a drop in the amplitude for 
the oscillation probabilities 1 — P{v^ v^) and P{v^ i/r), see Figs. I (b) and (d), 
respectively. Hence, the oscillation channel — > relevant for the atmospheric 
neutrino anomaly is influenced indirectly by matter effects through the enhancement 

of 013." 

In similarity to the two neutrino scheme, we found no sizable differences for the 
mantle-core-mantle model and the realistic Earth matter density profile. 

We also investigated the nadir angle dependence in the three neutrino case. 
The results are shown in Fig. ^ for different nadir angles ry. Note that for increasing 
nadir angle, the resonance for the core shrinks and the mantle starts to dominate the 

"In two neutrino analyses, the 1/^ Ur oscillation is usually expected not to be influenced by 
matter effects as there is no coupling to matter via the charged current interactions, which are 
present only for the electron flavor. 
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Fig. 3. Different transition probabilities as functions of the neutrino energy Ev for the realistic 
Earth matter density profile, (a) 1 — P{i'e —* ^e), (b) 1 — P(yij. — > I'lx), (c) P{i>a — > ^ij.,t), and (d) 




Fig. 4. The probability 1 — P{v^ — ► i/g) as a function of the neutrino energy E^, for different nadir 
angles. 
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resonance spectrum. For angles exceeding i] ~ 34° neutrinos do no longer traverse 
the core. The resonance is then centered around the resonant energy Ei, ~ 8.5 GeV 
corresponding to the matter density of the mantle (see Fig. |^ for rj — 45°). For 
larger nadir angles also the mantle resonance fades away due to the fact that the 
overall traveling path length of the neutrinos in the mantle decreases to zero. 

To achieve a large and well defined effect of matter enhancement, we thus propose 
to further study nadir angles just above the threshold where the traveling length is 
maximal through the mantle of the Earth. 

4. Atmospheric Neutrino Simulation 

To describe the quality of the step function approximation in a more intuitive way, 
we simulated a typical atmospheric neutrino experiment (like the S upe r-Kamiokande 
experiment) using the Honda et al. atmospheric neutrino fluxes t2l. The result is 
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Fig. 5. The left plot shows the energy spectrum of the electron neutrino-like events in an atmo- 
spheric neutrino detector for vacuum oscillation (solid), the realistic Earth matter density profile 
(dashed), and the step function approximation (dotted). The right plot shows the error of the 
step function approximation for different Amij- 

shown in Fig. ^. The left plot shows the energy spectrum of the core-crossing neu- 
trino bin (cos?/ = 0.8 -f- 1.0) in a 45kt-detector during one year of running. The 
figure, which is based on sin^ 26*23 = 1, sin^ 26*13 = 0.10, and Am^ = 3.2 • 10"^ eV^, 
shows the cases: vacuum oscillation (solid), realistic matter profile (dashed), and 
step function approximation (dotted) . The enhancements compared to the vacuum 
case are due to resonant matter effects. The right plot tries to quantify the difference 
of the two matter profiles. It shows the "error" of the step function approximation 
in the calculation of the total number of electron neutrino-like events in the core- 
crossing bin. The error is defined as the ratio K^Vrmf — Nsfa) / [N-rm-'p — ^VAc)|, 
where TVrmp, -^sfa, and AVac are the total number of events in the core-crossing 
bin for the realistic matter profile (RMP), the step function approximation (SFA), 
and vacuum oscillation (VAC), respectively. The defined quantity describes the 
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error of the step function approximation compared to the absolute magnitude of 
the matter effects. The error made using the step function approximation is at the 
level of 5%. Note that the total matter effects in the electron neutrino-like event 
rate of such an experiment ar e itself small (~ 10%). However, future muon storage 
ring neutrino experiments Eil promise better access to matter effects With such 
experiments it will possibly be necessary to take into account details of the matter 
density profile of the Earth. 



5. Summary and Conclusions 

We have studied the quality of a mantle-core-mantle step function approximation 
for the realistic Earth matter density profile in matter enhanced neutrino oscilla- 
tions. We considered a two neutrino and a three neutrino framework and found that 
the results given by the model function give an excellent approximation to the re- 
sults obtained with the realistic profile. The three neutrino analysis we performed, 
allowed us to investigate effects on the z/^ — > Vr channel, which are important for 
the description of the atmospheric neutrino anomaly. We found a drop of the cor- 
responding oscillation probability in the resonance region. A nadir angle dependent 
analysis of the enhancement showed that a sizable and clear resonance can be ex- 
pected for nadir angles around 34°, where the neutrinos spend maximal traveling 
length in the mantle and do not enter the core. Finally, we tried to quantify the 
quality of the step function approximation on the level of rates within a simulation 
of an atmospheric neutrino experiment. The obtained error is at the level of 5%. 
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